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ABSTRACT 

After at least 6 years of quiescence, Anomalous X-ray Pulsar (AXP) 4U 0142+61 entered an active 
phase in 2006 March that lasted several months and included six X-ray bursts as well as many 
changes in the persistent X-ray emission. The bursts, the first seen from this AXP in >11 years 
of Rossi X-ray Timing Explorer monitoring, all occurred in the interval between 2006 April 6 and 
2007 February 7. The burst durations ranged from 8— 3xl0 3 s. The first five burst spectra are 
well modeled by blackbodies, with temperatures kT ^ 2 — 6 keV. However, the sixth burst had a 
complicated spectrum that is well characterized by a blackbody plus three emission features whose 
amplitude varied throughout the burst. The most prominent feature was at 14.0 keV. Upon entry 
into the active phase the pulsar showed a significant change in pulse morphology and a likely timing 
glitch. The glitch had a total frequency jump of (1.9=t0.4)xl0“ 7 Hz, which recovered with a decay 
time of 17+2 days by more than the initial jump, implying a net spin-down of the pulsar. We discuss 
these events in the context of the magnetar model. 

Subject headings: — stars: neutron — X-rays: stars — X-rays: bursts — pulsars: individual 
(4U 0142+61) 


1. INTRODUCTION 

It is now generally accepted that the class of objects 
referred to as “Anomalous X-ray Pulsars” (AXPs) are 
magnetars - young isolated neutron stars powered by 
their high magnetic fields (Thompson &; Duncan 1995; 
Thompson & Duncan 1996). AXPs are X-ray pulsars 
wfith periods in the range 2-12 s, and X-ray luminosi- 
ties (~ 10 33 — 10 35 erg s -1 ) that cannot be accounted 
for by their available spin- down energy. The magnetar 
model was first proposed to explain the dramatic behav- 
ior exhibited by an apparently different object class - 
the Soft Gamma Repeaters (SGRs). SGRs show per- 
sistent properties similar to AXPs, but they were first 
discovered by their enormous bursts of soft gamma rays 
(> 10 44 erg) and their much more frequent, shorter, and 
thus less energetic bursts of hard X-rays. To date, SGR- 
like X-ray bursts have been observed from five AXPs, 
thus solidifying the connection between the two source 
classes (Gavriil et al. 2002; Kaspi et al. 2003; Woods 
et al. 2005; Krimm et al. 2006). For reviews of magnetar 
candidates and AXPs see Woods &; Thompson (2006), 
Kaspi (2007), and Mereghetti (2008). 

Thus far, only the magnetar model can explain the 
bursts observed from SGRs and AXPs (Thompson & 
Duncan 1995). The internal magnetic field exerts stresses 
on the crust which can lead to large scale rearrangements 
of the external field, which we observe as giant flares. If 
the stress is more localized, then it can fracture the crust 
and displace the footpoints of the external magnetic field 
which results in short X-ray bursts. The highly twisted 
internal magnetic field also slowly twists up the exter- 
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nal field; magnetospheres of magnetars may therefore be 
globally twisted (Thompson et al. 2002). Reconnection 
in this twisted magnetosphere has also been proposed as 
an additional mechanism for the short bursts (Lyutikov 
2002 ). 

In addition to bursts, AXPs and SGRs exhibit pulsed 
and persistent flux variations on multiple different 
timescales. An hours-long increase in the pulsed flux 
has been seen to follow a burst in AXP IE 1048.1-5937 
(Gavriil et al. 2006), On longer timescales, AXPs can 
exhibit abrupt increases in flux which decay on several- 
week timescales. These occur in conjunction with bursts 
and have been suggested as being due to thermal ra- 
diation from the stellar surface after the deposition of 
heat from bursts. Such flux enhancements have been 
observed in SGRs (see Woods et al. 2001, for exam- 
ple). The flux enhancement of AXP IE 2259+586 dur- 
ing its 2002 outburst can also interpreted as burst af- 
terglow (Woods et al. 2004), however, a magnetospheric 
interpretation is also plausible (Zhu et al. 2008). AXP 
IE 1048.1-5937 exhibited three unusual flux ‘flares. 5 In 
the first two, the pulsed flux rose on week-long timescales 
and subsequently decayed back on time scales of months 
(Gavriil & Kaspi 2004; Tam et al. 2008). These vari- 
ations have been tentatively attributed to twists im- 
planted in the external magnetosphere from stresses on 
the crust imposed by the internal magnetic field. AXPs 
XTE J1810-197, IE 1547.0-5408, and the AXP candi- 
date AX J1 845-0258 have also exhibited large flux vari- 
ations (Ibrahim et al. 2004; Halpern et al. 2008; Got- 
thelf & Vasisht 1998; Tam et al. 2006), however it is 
not clear whether these were of the abrupt rise type as 
in IE 2259+586 or the slow r -rise type as in IE 1048.1- 
5937. CXO J164710. 2— 455216 showed a clear abrupt 
rise (Muno et al. 2007). 1RXS J1708-4009 has been ar- 
gued to have flux variations associated with timing events 
(Israel et al. 2007b). Finally, AXP 4U 0142+61 has ex- 
hibited the longest timescale flux variations thus far, in 
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Fig. 1. - Epochs of the RXTE observations of 4U 0142+61 
analyzed in this paper (see § 3.2). Each observation is represented 
by a vertical line. The dashed line marks the entry of the source 
into the active phase. An “x” marks the first observation where the 
pulse profile was significantly different from the long-term average, 
“bl”, “b2”, and “b3” mark the observations containing bursts. To 
the left of the dotted line, the letters refer to groups of observations 
having similar total integration times. To the right, the letters refer 
to groups of observations having similar pulse profiles. 

which the pulsed flux increased by 29+8% over a period 
of 2.6 years (Dib et al. 2007; Gonzalez et al. 2009). 

4U 0142+61 is an 8.7-s AXP. It has a period deriva- 
tive of P ~ 0.2 x 10" 11 , implying a surface dipole mag- 
netic field of 1.3 x 10 14 G. 4U 0142+61 was monitored 
by RXTE in 1997 and from 2000 to 2007. Gavriil & 
Kaspi (2002) showed that 4U 0142+61 generally rotates 
with high stability. Morii, Kawai, k. Shibazaki (2005) 
reported a possible timing glitch in 1999 on the basis 
of an Advanced Satellite for Cosmology and Astrophysics 
(. ASCA ) observation in which the value of the frequency 
was marginally discrepant with the frequency as reported 
by Gavriil k Kaspi (2002). Dib et al. (2007) showed that 
the glitch may have occurred but is not required by the 
existing data. Dib et al. (2007) also reported on the evo- 
lution of the properties of 4U 0142+61 from 2000 March 
to 2006 April. In particular they reported stable timing, 
and an evolution of the pulse profile in 2—4 keV where 
the dip between the two peak was rising between 2000 
and 2006. They also reported a 29+8% increase in the 
pulsed flux between 2002 May and 2004 December. As of 
2006 March, in the published flux history of this source, 
there had been no reports of any X-ray activity such as 
bursts or flares, as described above for other AXPs. 

Here we report on the first detection of bursts from 
AXP 4U 0142+61, making this the sixth AXP for which 
this phenomenon has been observed. We also report that 
the source appears to have entered an active phase in 
2006 March in which almost every aspect of the emission 
changed. Our observations are described in Section 2. 
Our burst, pulsed morphology, pulse phase, pulsed flux, 
and timing analysis are presented, respectively, in Sec- 
tions 3.1, 3.2, 3.3, 3.4, and 3.5. In Section 4, we discuss 
the possible origins of this behavior and the implications 
for the magnet ar model. 

2. OBSERVATIONS 

4U 0142+61 has been monitored with Rossi X-ray 
Timing Explorer (RXTE) in 1997 and from 2000 to 2007. 
It has been monitored bi-monthly since 2005 March, with 
a typical observation length of ~5 ks. On 2006 March 
23, the source entered an active phase where many pulsed 
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TABLE 1 

Observations containing bursts. 


Observation 

ID 

MJD 

Date 

Number 
of Bursts 

Number of 
Active PCUs 

92006-05-03-00 

53831 

04/06/2006 

1 

3* 

92006-05-09-01 

53911 

06/25/2006 

4 

3 

92006-05-25-00 

54138 

02/07/2007 

1 

2 


One PCU switched ON and another one OFF partway through 
the observation, but the total number of active PCUs stayed con- 


stant throughout the observation. 

flux, spectral, and pulse profile changes were observed. 
We detected 6 bursts in three observations after the en- 
try into the active phase. After each burst detection, 
several RXTE target-of-opportunity (ToO) observations 
were made in addition to the regular monitoring. 

Here, we present a detailed analysis of the three obser- 
vations containing bursts (see Table 1). We also present 
an analysis of 91 other observations spanning the 2005 
March 21 to 2007 May 15 time period (MJD 53481 to 
MJD 54235, Observation IDs 91070-05-04-00 to 92006- 
05-32-00). 32 of these observations were prior to the 
entry into the active phase, and the remaining 59 ob- 
servations were after. This long-term analysis was per- 
formed in order to track the evolving pulsed flux and 
timing properties of the source. 

Figure 1 shows a timeline of the 94 analyzed observa- 
tions. Note that two segments of any observation that 
was split counted as separate observations if the seg- 
ments were given different observation IDs. The ranges 
of epochs with an increased density of observations con- 
tain the ToO observations. Prior to the active phase, 
the groups of observations referred to by capital letters 
in the Figure have similar total integration times. The 
groups of observations in the active phase have similar 
pulse profiles. These groups will be referred to in Sec- 
tions 3.2 and 3.5. 

3. ANALYSIS AND RESULTS 

All data presented here are from the Proportional 
Counter Array (PCA; Jahoda et al. 2006) aboard 
RXTE. The PCA is made up of five identical and in- 
dependent proportional counter units (PCUs). Each 
PCU is a Xenon/methane proportional counter with 
a propane veto layer. The data were collected in ei- 
ther GoodXenonwithPropane or GoodXenon mode which 
record photon arrival times with ~l-/zs resolution and 
bins them with 256 spectral channels in the ~2-60 keV 
band. 


3.1. Burst Analysis 

For each monitoring observation of 4U 0142+61, using 
software that can handle the raw telemetry data, we gen- 
erated 31. 25- ms lightcurves using all Xenon layers and 
events in the 2-20 keV band. These lightcurves were 
searched for bursts using our automated burst search al- 
gorithm introduced in Gavriil et al. (2002) and discussed 
further in Gavriil et al. (2004). In an observation on 2006 
April 6, we detected a significant burst, and four more 
bursts were detected in a single observation on 2006 June 
25. The sixth and most energetic burst w+is detected on 
2007 February 7. The bursts were significant in each ac- 
tive PCU. See Table 1 for the number of active PCUs in 
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each burst observation, as well as for the burst observa- 
tion epochs. 

To analyze these bursts we created event lists in FITS 4 5 
format using the standard FTQQLS 0 . For consistency with 
previous analyses of SGR/AXP bursts we extracted 
events in the 2-60 keV band. These events were re- 
duced to the solar system barycenter using the position 
found by Patel et al. (2003) for the source. The burst 
lightcurves are displayed in Figure 2 (LEFT). Our tech- 
niques for characterizing the temporal and spectral prop- 
erties of bursts were discussed in detail in Gavriil et al. 
(2004) but we repeat them here as some required modi- 
fication in order address the distinct properties of these 
bursts. 

Before measuring any burst parameters we deter- 
mined the instrumental background using the FT00L 
pcabackest. We extracted a background model 
lightcurve using the appropriate energy band and num- 
ber of PCUs. pcabackest determines the background 
count rate only on 16-s time intervals, so we interpolated 
these values by fitting a polynomial of order 6 to the en- 
tire observation, which yielded a good fit for each data 
set. 


3.1.1. Burst Temporal and Energetic Properties 

The burst peak time (£ p ), rise time (t r ) and peak flux 
Uv) were determined using the methods described in 
Gavriil et al. (2004). Usually, to measure the fluence 
for SGR and AXP bursts, we subtract the instrumental 
background for the lightcurve, integrate the light curve 
and fit the result to a step function with a linear term 
whose slope is the “local” background rate. The burst 
fluence in this case is the height of the step function. 
Although this technique worked well for the first burst, 
which was a short, isolated event, it was not appropriate 
for bursts 2, 3, and 4 because they had overlapping tails, 
and bursts 5 and 6 had tails that extended beyond the 
end of the observation. Thus, we opted to fit the burst 
flux time series, /(t), with exponential rises and decays: 
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The burst peak time (t p ), the burst rise time (£ r ), the 
peak flux (/ p ), the instrumental background model (b(t)) 
determined using pcabackest and the “local” back- 
ground (B) were held fixed at the values determined by 
the methods outlined above. We fit for the flux of the 
long-term tail component of the burst, and for the short 
(t/) and long-term (r t ) decay times of the tail of the 
burst. We then integrated our events to obtain fluence 
lightcurves. We fit the integral of Eq. 1 to our integrated 
events. The fits were performed over the entire span of 
the observations. Figure 2 (LEFT) shows our observed 
and model burst lightcurves; Figure 2 (RIGHT) shows 
our observed and model burst fluence time series. 

As is done for 7-ray bursts and SGR and AXP bursts, 
we characterized the burst duration by T90, the time from 
when 5% to 95% of the total burst counts have been 
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collected. To determine the X90 duration we integrated 
our burst model and numerically determined the 5% and 
95% time crossings. All burst temporal parameters are 
presented in Table 2. 

3.1.2. Burst Spectral Properties 

Burst spectra were extracted using all the counts 
within their Tgo interval. Background intervals were 
extracted from long, hand-selected, featureless intervals 
prior to the bursts. Response matrices were created us- 
ing the FTOOL pcarsp. The burst spectra were grouped 
such that there were at least 15 counts per bin after back- 
ground subtraction. Burst spectra, background spectra, 
and response matrices were then read into the spectral 
fitting package XSPEC 6 vl2.3.1. The spectra were fit to 
photoelectrically absorbed blackbodies using the column 
density found by Durant &, van Kerkwijk (2006). Only 
bins in the 2-30 ke-V band were included in the fits. The 
blackbody model provided an adequate fit for bursts 1 
through 5; model parameters are presented in Table 2. 

Extracted over its very long T 90 interval (see Table 2), 
the spectrum of Burst 6 was well fit by a simple black- 
body, however, there were hints of spectral features. On 
much shorter timescales, the spectrum of burst 6 was 
not well modeled by any simple continuum model and 
the presence of spectral features is clear (see Fig. 3). To 
study the temporal evolution of these features we ex- 
tracted 5-s long spectra in steps of 1 s from the peak 
of the burst, and repeated the spectral fitting procedure 
outlined above. Fig. 4 displays the spectral evolution 
of the features. The middle panel is a surface plot of 
the change in \ 2 after subtracting the continuum model 
(blackbody). A vertical slice in the middle panel corre- 
sponds to a spectrum of the burst extracted over a 5-s 
interval, and these spectra were extracted every 1 s . The 
features showed clear temporal variability and were most 
prominent near the onset of the burst (see Fig. 4). 

We looked closely at the first 5 s of the burst (see 
Fig. 3), and determined that a simple continuum model 
with less than three spectral features was not a statis- 
tically good representation of the data. A simple pho- 
toelectrically absorbed blackbody yielded a reduced \ 2 
of 2.54 for 41 degrees of freedom (DoF). Adding one or 
two Gaussian emission lines improved the fit but not suf- 
ficiently, yielding reduced x 2 s of 1.89 for 38 DoF, and 
1.67 for 35 DoF, respectively. Only a photoelectrically 
absorbed blackbody plus three Gaussian emission lines 
provided an adequate fit, with reduced x 2= =L21 for 32 
DoF. With this best- fit model, we find lines at energies 
8.8+J, 14.0+:! and 21.8 til keV - Our results from fit- 
ting different models to the first 5 s of the burst spec- 
trum are summarized in Figure 3, and the parameters 
returned from our best- fit model are listed in Table 3. 
We note that a power-law with index F~0.9lo!2 "together 
with three Gaussian emission lines also yielded a good 
fit, with reduced x 2 =1.20 for 32 DoF. In addition, us- 
ing a power-law instead of a blackbody as the continuum 
model returned values for the Gaussian emission line pa- 
rameters that were consistent with those listed in Table 3 
at the l-<7 level. 

From the surface plot in Fig. 4 (middle panel), we see 
that the 21.0 keV is only detectable at the 1-cr level at 

6 http://xspec.gsfc.nasa.gov 
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the start of the burst. The 14.0 keV is highly significant 
at the start of the burst, and remains detectable at the 
> 1-cr level for ~130 s. The 8.8 keV feature is the only 
one that remains detectable at the > l-o level beyond 
130 s from the burst onset. 

3.2. Pulse Profile Changes 

Many AXP outbursts are accompanied by significant 
pulse profile changes (e.g. Kaspi et al. 2003). To search 
for these in 4U 0142+61, for each observation, we gen- 
erated 64-bin pulse profiles using the method described 
in Dib et al. (2007). We then aligned the profiles with 
a high signal-to-noise template using a cross-correlation 
procedure in the Fourier domain. Then, for each group of 
observations in the active phase, we summed the aligned 
profiles, extracted the DC component from the summed 
profile, and scaled the resulting profile so that the value 
of the highest bin is unity and the lowest bin is zero. 

The resulting pulse profiles are shown in Figure 5 in 
chronological order. The panels in the left-most column 
are from Dib et al. (2007) and show pulse profiles in the 
6 years prior to the entry into the active phase. Notice 
the slow change in the height of the dip between the 
peaks. The panels in columns 2 and 3, marked with the 
corresponding letter in the top right corner, show pulse 
profiles for each of the data groups in the active phase 
that were marked with a letter in Figure 1. The 4 plots 
in bold correspond to the observations marked with “x” , 
“bl”, “b2”, and “b3” in Figure 1. 

The pulse profile evolution can be described as fol- 
lows: prior to the active phase, features in the double 
peaked pulse profile were evolving on a timescale of sev- 
eral years (see panels C6 to CIO corresponding to RXTE 
Cycles 6 to 10). Then, in the first observation of the ac- 
tive phase (panel x), the pulse profile became suddenly 
triple-peaked. It was also triple-peaked in the following 
observation (panel bl), in which a burst occurred. It re- 
mained multi-peaked for the following two groups of ob- 
servations (panels G and H). Then, an observation with 
multiple bursts occurred (panel b2) . In that observation, 
also having a triple-peaked pulse profile, the middle peak 
was taller than the other two. Following this observation, 
the pulse profile seemed to be slowly recovering back to 
its double-peaked long-term shape (panels I, J, K, and 
L). Another burst observation interrupted this evolution 
seven months later (panel b3). In that observation, a 
large burst was detected. In the pulse profile of that ob- 
servation, the left-most peak was significantly taller than 
usual. The event that caused this change had apparently 
no effect on the following observation which occurred 2 
days later: the profile went back to being double-peaked 
(panels M and N). To summarize, the pulse profile be- 
came multi-peaked at the beginning of the active phase. 
Following the second observation with bursts, the profile 
started to recover to its double-peaked shape. The evo- 
lution was only temporarily interrupted for the duration 
of the third observation with bursts. The behavior of the 
pulse profile in the 2—4 and 4—10 keV bands was similar. 

Note that from Figure 5 alone we can compare the 
sizes of the peaks to each other, but we cannot track the 
evolution of the heights of each peak separately. In order 
to do that, we must scale the pulse profile of each group 
of observations by the average pulsed flux of that group. 
This analysis is presented in Section 3.4.3. 


We also performed an analysis of the Fourier compo- 
nents of the pulse profiles. The results are shown in Fig- 
ure 6. The variations in the power of the fundamental 
Fourier component are shown in panel 1, and that of the 
second harmonic in panel 2. Note how the amplitude of 
the fundamental varied monotonically prior to, but not 
during, the active phase. Also note how the power in the 
second harmonic was already back to its pre-active-phase 
level before the last burst occurred. 


3.3. Burst Rotational Phases 

An important factor in understanding the origin of the 
bursts is the rotational phase at which they occur. The 
phases of the bursts are shown in Figure 7. For each 
burst observation we created at 31.25-ms time resolution 
lightcurve and folded it using our timing solution (see 
§ 3.5). We then phase-aligned these folded profiles by 
cross-correlating them with the long-term pulse profile 
template. Our phase-aligned folded profiles are shown 
in Figure 7 (histograms in the last row). Superposed 
on each folded profile are two curves. The top curve is 
made of the 5 Fourier harmonics that best fit the his- 
togram. The bottom curve is made of the best-fit 5 har- 
monics after the removal of the 4 seconds centered on 
each burst. Note how the two curves in a given plot 
are similar, demonstrating that the additional peaks in 
the profiles are not due to the burst. The first three 
bursts occurred near the middle “new” peak in the pro- 
file. Burst 5 occurred near the old tall peak of the profile 
(see Fig. 5). The phase of burst 6 corresponds again to 
a new peak in the profile, this time located where the 
previous small peak used to be. The coincidences of sev- 
eral bursts with new, transient profile features that are 
present even when the actual burst data are removed 
are suggestive of lower-level transient emission from the 
same physical location, with the burst being the extreme 
of this emission’s luminosity distribution. 


3.4. Pulsed Flux Analysis 
3.4.1. Short-Term Pulsed Flux Analysis 

Previous AXP bursts are often accompanied by short- 
term pulsed flux enhancements (e.g. Gavriil et al. 2006). 
To search for these, for each of the three observations 
containing bursts, we made two barycentered time series 
in count rate per PCU, one for the 2—4 keV band and 
the other for 4—20 keV. We included only photons de- 
tected by PCUs that were on for the entire duration of 
the observation. The time resolution was 31.25 ms. We 
removed the 4 s centered on each burst from each time 
series. Then, we broke each time series into segments 
of length ~500 s. For each segment, we calculated the 
pulsed flux using two different methods. 

First, we calculated the RMS pulsed flux using 


RMS 


- , 2 £«« : 
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where is the even Fourier component defined as 
a k = jjYliLiPi cos (Z'Kki/N), cr ak is the uncertainty of 
afc, bk is the A; th odd Fourier component defined as bk 

= ^J2iLiPi s * n (27rki/N), a t, k is the uncertainty of bk, i 
refers to the phase bin, N is the total number of phase 
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TABLE 2 

4U 0142+61 Burst Temporal and Spectral Properties 


Parameter 3, 

April 2006 


June 2006 


February 2007 

Burst 1 

Hurst 2 

Burst 3 

Burst 4 

Burst 5 

Burst 6 

Burst day (MJD UTC) 
Burst start time (UT) 
Rise Time b 

t 90 c 

Phase d 
Tqq Fluence e 
Tqo Fluence f 
Peak Flux g 
Peak Flux h 
kT' 1 

53831 

07:09.55.544(7) 

7 + 4 
8.0±^ 
0.503(7) 

65+6 
16+3 
299+59 
75+19 
6 2 +4 - 4 

+£-=1,9 

53911 

• 01:15:54.555(11) 

+3 

562.2+1!’ 

0.487(4) 

1847+4 

229+11 

367+65 

46+8 

2 R7+ 0 - 10 
z ' D Enng 

53911 

01:15:55.119(43) 

43^22 

6 44° 15 
D - 44 -0.15 

0.54(2) 
125+9 
26+3 
250+54 
53+12 
4 SO 4 ” 0-38 

53911 

01:16:9.216(4) 

4+ 

H.62++ 

0.173(4) 

147+8 

38+3 

228+52 

58+14 

5 16+ 0 - 44 
++P-0.34 

53911 

01:20:0.131(3) 

3± 3 

92.54^0 g5 

0.749(4) 

342+9 

43+4 

245+54 

31+7 

9 c;9+°- 2 1 
£+-£=0.13 

54138 

10:04:43.264(27) 

27± 4 

3472.9^1 4 
0.358(4) 
14377+41 
1455+53 
551+97 
253+22 
9 1 O+0-98 

z :l z - 0.01 


a All quoted errors represent 1 -a uncertainties. 
b Rise time (ms). 

c T90 duration (s), defined as the time it took to collect 90% of the total burst fluence. 
d Phase (cycles), phases relative to the template shown in Figure 7. 
e 2-60 keV fluence in the Tqq duration in units of counts PCU” 1 . 
f 2-60 keV fluence in the Tqo duration in units of 10" 10 erg cm” 2 . 
g 2-60 keV peak flux for 64 ms in units of counts s” 1 PCU” 1 . 
b 2-60 keV peak flux for 64 ms in units of 10~ 10 erg s” 1 cm” 2 . 

Blackbody temperature (keV). For burst 6, the temperature is that of the blackbody component after accounting for the apparent 
spectral features (see §3.1.2). 
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Pig. 2. LEFT - 2-60 keV burst lightcurves binned with 31.25-ms time resolution as observed by RXTE. The red curves are the best-fit 
exponential rise and decay models (Eq. 1). RIGHT - 2-60 keV integrated burst lightcurves. For each burst lightcurve we subtracted the 
instrumental background, the background due to other sources in the FOV, and the contribution of the tails of neighboring bursts before 
integrating. Notice that some bursts have very long tails. The red curve is the integral of the model in the left panel. The blue curve is 
a step function whose height corresponds to the location of the total fluence as determined by extrapolating the best-fit integrated burst 
model (red curve) to infinity. 
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Fig. 3. — Panel 1A: The spectrum of burst 6 extracted over 
a 5 s interval starting from 1 s before the burst peak. The red 
solid lines represents the best-fit model- a blackbody (blue curve) 
plus 3 Gaussian emission lines (green curves), all photoelectrically 
absorbed. The fit had a reduced x 2 of 1.2 for 32 degrees of freedom 
(DoF). Panel IB: Residuals, expressed in terms of the change in % 2 , 
after subtracting the best-fit model (blackbody 4 - lines) shown in 
Panel 1A. Panels 2 A and 2 B : Same, except here the best-fit model 
(red curve) consists of a blackbody (blue curve) plus 2 Gaussian 
emission lines (green curves), all photoelectrically absorbed. The 
fit had a reduced \ 2 of 1.67 for 35 DoF. Panels 3 A and 3 B: Same, 
except here the best-fit model (red curve) consists of a blackbody 
(blue curve) plus a single Gaussian emission lines (green curve), 
both photoelectrically absorbed. The fit had a reduced x 2 of 1.89 
for 38 DoF. Panel 4A and 4B: Same, except here the best-fit model 
consists solely of a photoelectrically blackbody (red curve). The 
fit had a reduced x 2 of 2.54 for 41 DoF. 

bins, pi is the count rate in the phase bin of the 
pulse profile, and n is the maximum number of Fourier 
harmonics to be taken into account; here n— 5 for consis- 
tency with Dib et al. (2007) and Gonzalez et al. (2009). 

While least sensitive to noise compared to other pulsed 
flux measurement methods, the RMS method returns 
a pulsed flux number that is affected by pulse profile 
changes. To confirm our pulsed flux results, we also used 
an area- based estimator to calculate the pulsed flux: 



*= i 


( 3 ) 


where p m i n is the average count rate in the off-pulse phase 
bins of the profile, determined by cross-correlating with 
a high signal-to-noise template, and calculated in the 
Fourier domain after truncating the Fourier series to 5 


TABLE 3 

Burst 6 Spectral Fit 


Parameter 

Value a 

Reduced x 2 (Degrees of Freedom) 
2-60 keV Flux b (10 ~ 9 erg s~ 1 cm“ 2 ) 
2-20 keV Flux b (10 ~ 9 erg s^cm- 2 ) 

1.21(32) 

o qo + 0.48 

O A' 7 — 0.12 
2 , 67 - 0.06 

Blackbody Component 

Temperature, kT (keV) 

Radius 0 (km) 2.3 

2-60 keV Flux b (10 “ 9 erg s~ 1 cm“ 2 ) 

2-20 keV Flux b (10 “ 9 erg s~ 1 cm -2 ) 

1 S +°- 4 
i - S - 0.3 

iSMofc) 

0.02 

0 ^q + 0 - 02 
u -°^-o.oi 

Gaussian Emission Line 1 

Energy (keV) 

Width, u (keV) 

2-60 keV Flux b (10 ™ 9 erg s“ 1 cm~ 2 ) 
2-20 keV Flux b ( 10“ 9 erg s - 1 cm" 2 ) 

0 q + 0.2 
— 0.2 

U ' y - 0.3 

0 250" 1 " 0 ' 031 
u - zou -o . 010 

0 250”f 9,9 ^ 
u.^ou _ 0 006 

Gaussian Emission Line 2 

Energy (keV) 

Width, cr (keV) 

2-60 keV Flux b (10 ~ 9 erg s - 1 cm“ 2 ) 
2-20 keV Flux b (10 “ 9 erg s~ 1 cm“ 2 ) 

14 fC 0 ' 4 
9 0 + O.4 

1 c ' 7 + 0.19 
J -- J * - 0.06 
J-*oo _ 0 04 

Gaussian Emission Line 3 

Energy (keV) 

Width, cr (keV) 

2-60 keV Flux b (10 ~ 9 erg s“ 1 cm~ 2 ) 
2-20 keV Flux b (10 -9 erg s _ 1 cm -2 ) 

21 . 8 i^ 

0 n+1-7 
•J-U — 0 .g 

1 * 615 - 0.16 
n qchFO - 016 

0.360_q 008 


a All errors were extracted at the 1 cr level. 

b All quoted fluxes are unabsorbed. In the fit, the column density 
was held fixed at N H ~ 0.64x10 22 cm “ 2 (Durant & van Kerkwijk 
2006). 

c The blackbody radius was calculated assuming a distance to the 
source of 2.5 kpc (Hulleman et al. 2004). 

harmonics. The results are shown for Frms in Figure 8 
for F rms (FArca gives consistent results). Note the sig- 
nificant increase in the 4-20 keV pulsed flux in the 2006 
June observation following the cluster of bursts. This 
increase is not present in 2—4 keV. Also note the signifi- 
cant rise and subsequent decay of the pulsed flux follow- 
ing the large 2007 February burst. The pulsed flux was 
sufficiently enhanced in these two observations that one 
can see individual pulsations by eye in Figure 7 in the 
two panels marked with a star, containing the raw burst 
lightcurves with 31.25-ms time resolution. 

3.4.2. Long-Term Pulsed Flux Analysis 

For each of the 94 analyzed observations, we created a 
pulse profile (in units of count rate per PCU) using the 
same procedure as in Section 3.2. Then we calculated 
the pulsed flux for each observation using Equations 2 
and 3. Data from PCU 0 were omitted because the long- 
term trend in the pulsed counts is not the same as that 
in the remaining PCUs, presumably due to the loss of its 
propane layer. Data from PCU 1 were omitted after the 
loss of its propane layer as well, on MJD 54094. 

We extracted pulsed fluxes in the 2—4 and 4—20 keV 
bands using both the RMS and area pulsed flux method 
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Fig. 4. — Spectral evolution of burst 6. Top: Burst lightcurve binned with 1-s time resolution and extracted over the same energy 
band used for spectral fitting, 2-30 keV. Middle: Surface plot of the difference between the burst spectrum and a simple continuum model 
(blackbody) in terms of change in x 2 ■ Each vertical slice in the surface plot corresponds to a spectrum that was calculated over a 5-s long 
window. This window was translated across the burst in 1-s long steps. For each spectrum, a best-fit blackbody model was subtracted. The 
color wedge on the right maps the colors to the resulting change in \ 2 • The contours indicate the regions where the spectrum differed by 
the model by 1-a (68.3% confidence, blue contour), 2-cr (95.4% condfidence, green contour) and 3-a (99.7%, yellow contour), respectively. 
Notice how the feature at 21.8 keV (see Table 3) is detected at the 1 -<t level only at the start of the burst. The feature at 14.0 keV is highly 
significant at the start of the bursts, and remains detectable at the >l-cr level for ~130 s. Notice, that beyond ^130 s the only feature 
listed in Table 3 that remains detectable at the >l-a level is the feature at 8.8 keV. Bottom: The insets are examples of individual 5-s 
long spectra used to generate the surface plot. These spectra were extracted at 0, 50 and 140 s from the burst onset, respectively. The red 
curves indicate the corresponding best-fit blackbody model. 


because of the numerous pulse profile changes around the 
times of the bursts. There is no evidence of a long-term 
change in the pulsed flux associated with the bursts ex- 
cept, possibly, for a hint of an increase in the 2—4 keV 
band roughly midway between the second and third ob- 
servations containing bursts. We also note that the 
4—20 keV pulsed flux of two of the observations con- 
taining bursts are significantly larger than the long-term 
average. Removing the bursts from these observations 
does not change this result. 

We performed the same analysis for individual obser- 
vations in 2—10 keV for an extended period of time. This 
is shown in panels 2 and 3 of Figure 9. This plot is an 
update to Figure 10 of Dib et al. (2007). 

In order to verify that trends seen in panels 2 and 3 
of Figure 9 are not an artifact of the evolution of the 
response of the detector, we calculated the pulsed flux 
in erg s" 1 cm” 2 using a method that took the evolution 
of the response into account. The method is described 
in Dib et al. (2007) and takes into account spectral fits 


obtained from imaging data (in this case XMM data, 
see Gonzalez et al. 2009) to convert counts to energy for 
each combined set of observations. For this analysis, we 
used data from all PCUs; however, data from PCU 1 
were excluded after the loss its propane layer. Data from 
PCU 0 were included because the response matrices used 
took into account the loss of its propane layer. 

The results are shown in panel 1 of Figure 9. The 
first 6 points, corresponding to RXTE Cycles 6—10, are 
from Dib et al. (2007). The second to last point is ob- 
tained by combining all observations that occurred be- 
tween bursts 1 and 2, but omitting observations contain- 
ing bursts. The observations that we included took place 
during the exponential recovery of the possible glitch (see 
Section 3.5). The last point in panel 1 was obtained by 
combining the observations that occurred after burst 2, 
but again omitting those containing bursts. The obser- 
vations we included took place after the end of the expo- 
nential recovery of the possible glitch (see Section 3.5). 

In the first of the two data points in the active phase, 
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Fig. 5. Pulse profile evolution of 4U 0142+61. Left-most column: normalized average 2-10 keV pulse profiles in the 6 years prior to 
the entry into the active phase in chronological order (panels from Dib et al. 2007). Middle and right-most columns: normalized, averaged 
2—10 keV pulse profiles for each of the data groups marked by a letter in Figure 1, after the entry into the active phase. The 4 plots with 
bold labels correspond to the observations marked with “x” , “bl”, L£ b2” , and “b3” in Figure 1. 


the pulsed flux in erg s -1 cm~ 2 is 18±8% larger than 
in the pre-active phase. This is consistent with the in- 
crease reported in Gonzalez et al. (2009) in the same 
energy range. A hint of this increase can be seen in 
panel 2 although it appears less significant. This discrep- 
ancy can be accounted for by the fact that the spectrum 
changed: Gonzalez et al. (2009) reported a temporary in- 
crease in the spectral hardness in an XMM observation 
of 4U 0142+61 immediately following the bursts. 

3.4.3. Combined Pulse Morphology and Pulsed Flux Analysis 

In Section 3.2, we calculated the Fourier components 
of the aligned average pulse profiles. This gave us the 
relative amplitude of the pulse profile harmonics in each 
group of observations marked with a letter in Figure 1. 
In Section 3.4.2, we calculated the pulsed flux for every 
observation. Here, we compute a weighted average of 
the pulsed flux for each group of observations using the 
flux points calculated in Section 3.4.2. We then recon- 
struct the profiles for each of the groups from the first 


five Fourier components (not including the DC), scale 
them by the average RMS pulsed flux for that group, 
and add the necessary offset for the lowest point on each 
curve to be zero. This means that the resulting scaled 
profiles return the correct pulsed flux. The advantage of 
this analysis is that we can now trace the evolution of 
each of the peaks independently. 

The results are presented in Figure 10. In panel (a), we 
show the average scaled pulse profile per observation for 
observations in the year preceding the active phase. The 
profile is double-peaked. In panel (b), we show in bold 
the scaled profile for the observation containing burst 1 . 
We also show the scaled profiles in each of the groups 
of observations that followed the burst. The profiles are 
triple-peaked. In panel (c), we show in bold the scaled 
profile for the observation containing bursts 2, 3, 4, and 5. 
We also show the scaled profiles in each of the groups of 
observations that followed the bursts. Note how the in- 
crease in the pulsed flux in the observation containing the 
bursts is not only a consequence of the appearance of the 



Activity from AXP 4U 0142+61 


9 


Year 


2002 2004 2006 



52000 52500 53000 53500 54000 


Modified Julian Day 


Fig. 6. - Top: evolution of the power in the fundamental Fourier component of the pulse profile of 4U 0142+61 in the 2—10 keV 
band. Bottom: Evolution of the power in the second harmonic. In both panels, the points with large horizontal error bars are from Dib 
et al. (2007). The remaining points are obtained from groups of observations after the entry into the active phase. The three bold points 
correspond to the observations with bursts. The three dashed lines correspond to burst epochs. 


new peak, but a result of the increase in size of all three with 64 phase bins, were cross-correlated in the Fourier 

peaks. Also note the evolution of the pulse profile back domain with a high signal-to-noise template created by 

to being double-peaked. In panel (d), we show in bold adding phase-aligned profiles from all observations. The 

the scaled profile for the observation containing burst 6. cross-correlation returned an average pulse time of ar- 

We also show the scaled profiles in each of the groups rival (TOA) for each observation corresponding to a fixed 

of observations that followed the burst. Note again how pulse phase. The pulse phase (f) at any time t can be ex- 

the pulsed flux increase is due to both peaks increasing pressed as a Taylor expansion, 

in size. 

3.5. Timing Analysis <p(t) ~ <po(to) + u 0 (t - t Q ) + -i/o(£ ~ *o) 2 

Many AXP outbursts and active phases are accompa- 1 

nied by interesting timing anomalies (e.g. Kaspi et al. — A>) + • • •? (4) 

2003; Woods et al. 2004; Israel et al. 2007a; Dib ° 

et al. 2009). Next we consider the timing behavior of where v = 1/P is the pulse frequency, v = du/dt, etc., 

4U 0142+61, which, as we show, also exhibits interesting an( j subscript ‘0” denotes a parameter evaluated at the 

evolution at the start of the active phase. ^ reference epoch t — to- To obtain ephemerides, the TOAs 

For all our R A TE observations of 4U 0142+61, photon were fitted to the above polynomial using the pulsar tim- 

arrival times at each epoch were adjusted to the solar j n g software package TEMPO 7 . 

system barycenter using the position obtained by Patel As explained above, our first, set of TOAs was obtained 

et al. (2003) from Chandra data. They were then binned by aligning the folded observations with a template pro- 

with 31.25-ms time resolution. In the timing analysis file using a cross-correlation procedure. In order to de- 

presented below, we included only the events in the en- termine to what extent this set of TOAs was affected by 
ergy range 2.5-9 keV, to maximize the signal-to-noise the pulse profile changes that took place in the active 

ratio of the pulse. phase, we generated two additional sets of TOAs (sets 2 

Each barycentric binned time series was epoch- folded an d 3^ 

using an ephemeris determined iteratively by maintain- 
ing phase coherence; see below. Resulting pulse profiles, 7 See http://www.atnf.csiro.au/research/pulsar/tempo. 
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Fig. 7. — Phase analysis of the 6 detected bursts. Each column corresponds to an observation in which bursts were detected. Top row: 
200-s segments of the time series containing the bursts in the 2-20 keV band. The time resolution is 1 s. Middle row: 4-s segments of 
the time series containing the bursts. The time resolution 31.25 ms. Bottom row: aligned folds of the burst observations shown below 
the scaled long-term average profile. The folded counts are presented as histograms. Superposed on each histogram are two curves. The 
top curve is made of the 5 Fourier harmonics that best fit the histogram. The bottom curve is made of the best-fit 5 harmonics after the 
removal of the 4 seconds centered on each burst. The arrows indicate the phases of each burst. 


For set 2, we made the assumption that the location of 
the off-pulse, determined by finding the minimum in the 
pulse profile after smoothing, is not affected by the pulse 
profile changes. We then generated TOAs by aligning 
the off-pulse regions of the folded observations with the 
off-pulse region of the template and extracting phase dif- 
ferences. The last set of TOAs (Set 3) was obtained after 
aligning the tallest peak of each smoothed pulse profile 
with that of the long-term template. 

Timing residuals for all three sets of TOAs are shown 
in Figure 11 using two different ephemerides. The linear 
ephemeris used in the top panel is that obtained by fitting 
the first set of TOAs in the regions indicated by the hori- 
zontal arrows, with a phase jump in between. The linear 
ephemeris used in the bottom panel is the ephemeris that 
best fit the first set of TOAs prior to the active phase. 
The same long-term trend is present in the first two sets 
of timing residuals, with more scatter in the second set. 
The difference in phase between each bold circle and the 
corresponding empty circle represents our uncertainty in 
determining a fiducial point on the pulsar. In the third 
set of residuals, the outlier points represent the observa- 
tions where the largest peak was no longer the right-most 
peak. Apart from the outliers, the same trend seen in the 
other two sets of residuals is present in the third set. 

From here on, we assume that the presence of the same 
trend in our three sets of residuals is an indication that 
the TOAs in the first set were not significantly affected 
by the pulse profile changes. We therefore have used the 


first set of TOAs in the remainder of this Section. 

Using the above assumption, the results of the tim- 
ing analysis are shown in Figure 12. The results of a 
segmented timing analysis are in panels 1 and 2. The 
results of a long-term timing analysis are in panels 3, 4, 
5, and 6 of the same Figure. 

For the segmented timing analysis, we used the 70 
TOAs that had the smallest uncertainties and omitted 
the TOAs of the 3 burst observations as well as the TOA 
of the first observation in the active phase. We then 
divided the data into segments of similar pulse profiles. 
For each data segment we found a linear ephemeris us- 
ing TEMPO 8 . We plotted the resulting ephemerides in 
panel 1 of Figure 12 with the uncertainties returned by 
TEMPO. The timing residuals are shown in panel 2. No- 
tice how the slope between bursts 1 and 2 is more neg- 
ative than the long-term average. 

We then fit a linear trend through all the pre-burst 
observations, and another linear ephemeris through the 
post-burst observations after the pulse profile had started 
to return to the double-peaked shape. We did not include 
the TOA corresponding to the last observation contain- 
ing bursts. The results are plotted in panel 3 of Figure 12 
with uncertainties, with the residuals in panel 4. Notice 
the difference in the slope in the two regions. In partic- 
ular, extrapolating the two ephemerides to a point be- 
tween their times of validity makes it seem like a sudden 

8 For the segmented timing analysis we ran TEMPO in mode 1. 
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Fig. 8. — RMS pulsed flux within the observations containing bursts. Each column corresponds to one observation. In each column we 
show, descending vertically, the 1-s resolution lightcurve with the bursts indicated, the 2-4 keV RMS pulsed flux, and the 4-20 keV RMS 
pulsed flux. The dotted line in each of the pulsed flux plots shows the average of the pulsed fluxes obtained after segmenting and analyzing 
the time series of the observation immediately prior to the one shown. 


spin down, i.e. an “anti-glitch” occurred. 

Finally, we included all the TOAs and did one global 
fit. In order to provide a good fit to the TOAs at the 
onset of the active phase, we had to assume that a glitch 
occurred on MJD 53809, with the glitch model consisting 
of a permanent change in m and z> and a frequency change 
A Vd that recovered exponentially on a timescale t^, i.e., 

V = 14) (t) + A 1/ + A v d e~ {t - t9)/Td + A J> (t - tg ), (5) 

where mo(£) is the frequency evolution pre-glitch given by 
m o(t) ~ i/(to) + v(t — t 0 ), Am is a instantaneous unrecov- 
ered frequency jump, A is the frequency increase that 
decays exponentially on a time scale t g is the glitch 
epoch determined by setting the phase jump to zero, and 
Az> is the post-glitch change in the long-term frequency 
derivative. The values of the fit parameters are listed in 
Table 4. 

From Table 4, we can see that the total sudden change 
in frequency which happened at the onset of the active 
phase (Avtot = Am + A m^) is positive, and Am^ decayed 
exponentially. At the end of the decay, the net effect 
of the recovered glitch was a negative Am. That is, the 
data suggest that after had decayed, there remained 
net spin- down relative to the undisturbed ephemeris from 
before the active phase. 

The fit parameter Am^, from which we conclude a sud- 
den spin-up, is controlled primarily by the first few TOAs 
in the active phase. It is therefore possible that it is af- 
fected by pulse profile changes. However, as shown in 
Figure 11, the TOAs in the active phase were not signif- 
icantly contaminated by pulse profile changes. Also, the 
segmented analysis shown in Figure 12 clearly indicates 
an initial spin-up. 

In contrast to Am^, the fit parameter Am, from which 
we conclude a net spin down, is primarily controlled by 


TABLE 4 

Spin and Glitch Parameters for 4U 0 142+6 l a 


Parameter 

Value 

MJD start 

53481.268 

MJD end 

54235.319 

TOAs 

93 b 

V (Hz) 

0.1150920955(12) 

V (lcr 14 Hz s- 1 ) 

-2.661(9) 

Epoch (MJD) 

53809.185840 

Glitch Epoch (MJD) 

53809.185840 

(Hz) 

— 1.27(17)xl0~ 8 

Az> (Hz s- 1 ) 

-3.1(1. 2)xl0~ 16 

Ai/ d (Hz) 

2.0(4) Xl0“ 7 

t d (days) 

17.0(1.7) 

RMS residual (phase) 

0.0168 


a Numbers in parentheses are TEMPO-reported la uncertainties. 
b A single TOA was omitted due to the very poor signal to noise 
ratio in the corresponding observation. 


the TOAs outside of the active phase, which were cer- 
tainly not affected by pulse profile changes (Figure 12, 
panel 3). 

We therefore conclude that the pulsar likely suffered a 
spin-up glitch near or at the start of the radiatively active 
phase, but that the glitch 1 over-recovered’ such that long- 
term, its net effect is a spin- down of the pulsar. 

4. DISCUSSION 

In this paper, we have described the timing, pulse pro- 
file, and pulsed flux behavior of 4U 0142+61, during its 
2006-2007 active phase, the first such episode yet stud- 
ied from this source. Specifically we have shown that in 
addition to a sudden departure from a slow, systematic 
evolution of the source’s pulsed flux and pulse profile, 
this AXP also suffered a significant timing event that is 
best described as a sudden spin- up glitch, followed by 
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Modified Julian Day 

Fig. 9. Long-term pulsed flux time series in 2—10 keV 
for 4U 0142+61. Panel 1: area pulsed flux in erg s -1 cm - 2 
for combined observations. The solid line marks the 29+8% in- 
crease reported in Dib et al. (2007). Panel 2: area pulsed flux in 
counts/s/PCU for individual observations. Panel 3: RMS pulsed 
flux in counts/s/PCU for individual observations. All panels: the 
dashed lines mark the burst epochs. The points marked with stars 
are the pulsed fluxes of the observations containing bursts. 

a large decay of the frequency jump such that the net 
effect was a slow-down with Ai *jv ~ —8 x 10" 8 . In- 
terestingly, the pulsed X-ray flux showed no significant 
change apart from immediately following bursts, in con- 
trast to other AXP radiative outbursts (e.g. Kaspi et al. 
2003; Tam et al. 2008; Ibrahim et al. 2004; Muno et al. 
2007). Further, we have reported on six SGR-like bursts 
from the source that occurred during this period, one of 
which was notable for its unusual spectrum, which was 
poorly fit by an continuum model. 

4.1. The SGR-Like Bursts from 4U 0142+61 

Five of the six bursts reported here for 4U 0142+61 had 
fast- rise-slow- decay profiles, with tails much longer than 
the rise times. These morphologies are similar to the 
class of AXP bursts labeled as “Type B” by Woods et al. 
(2005). Type B AXP bursts have also been character- 
ized by arriving preferentially at pulse maximum, and are 
seen more often in AXPs compared with SGR s (Woods 
et al. 2005; Gavriil et al. 2006). They were suggested as 
being due to a sudden rearrangement of magnetic field 
lines anchored in the crust following a crustal fracture 
(Thompson & Duncan 1995), as opposed to reconnection 
events in the upper magnetosphere (e.g. Lyutikov 2002), 
argued as more likely for the shorter, symmetric “Type 
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Fig. 10.- - Average pulse profile per observation in the 2—10 keV 
band in several groups of observations reconstructed from the first 
five Fourier components, and scaled to return the appropriate 
pulsed flux, (a) average scaled pulse profile for observations in 
RXTE Cycle 10, the last RXTE Cycle before the active phase. 
(b) Bold curve: scaled pulse profile of the observation containing 
burst 1. Thin curves: average scaled pulse profile for each of the 
groups of observations that followed burst 1 (groups G and H in 
Figure 1). (c) Bold curve: scaled pulse profile of the observation 
containing bursts 2, 3. 4, and 5. Thin curves: average scaled pulse 
profile for each of the groups of observations that followed the 
bursts (groups I, J, K, and L in Figure 1). (d) Bold curve: scaled 
pulse profile of the observation containing burst 6. Thin curves: 
average settled pulse profile for each of the groups of observations 
that followed burst 6 (groups M, and N from Figure 1). 
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Fig. 11. — Top: Timing residuals for all three sets of TOAs. Bold circles: timing residuals obtained using the first set of TOAs obtained 
by cross correlation in the Fourier domain. Empty circles: residuals obtained using the second set of TOAs obtained by aligning the 
off-pulse regions of the pulse profiles. Triangles: residuals obtained using the third set of TOAs obtained by aligning the tallest peak of 
the pulse profiles. The linear ephemeris used to produce all three sets of residuals is that which best fit the first set of TOAs in the regions 
indicated with the horizontal arrows (with an arbitrary phase jump in between). Bottom: Timing residuals for the same three sets of 
TOAs. The linear ephemeris used is that which best fit the first set of TOAs in the pre-active phase region. Both panels: the dotted line 
indicates the epoch of the first observation of the active phase. The dashed lines indicate the burst epochs. 


A” bursts that show no preference for pulse maximum. 
However, it is interesting that in the 4U 0142+61 bursts, 
in spite of their B-type morphology, no clear preference 
for arrival at or near pulse maximum was seen (see Table 
2). This may blur somewhat the distinction between the 
putative types, although the morphological distinction 
remains clear. 

It is also notable that the vast majority of bursts seen 
from three of the four AXPs from which bursts have yet 
been detected, IE 1048.1-5937, XTE J1810-197 and now 
4U 0142+61, had long tails following fast rises. The ex- 
ception todate is IE 2259+586 for which, during its 2002 
outburst, the minority (roughly one dozen out of 80) of 
the bursts seen were of this form, the majority being of 
Type A, similar to those classically seen in SGRs, and 
suggested by Woods et al. (2005) to be magnetospheric. 
A possible hint regarding the origin of the different types 
may lie in that IE 2259+586 was mid-outburst when its 
bursts were observed, whereas the other sources 7 bursts 
all occurred in the days/months following the commence- 
ment of an active period, presumably as the pulsar re- 
covered from a major event, rather than mid- event. 


4.1.1. Spectral Features 

The complicated spectrum of burst 6 is puzzling. The 
most significant feature at ~14 keV is particularly in- 
teresting. Emission features at similar energies were ob- 
served from two out of the three bursts from IE 1048.1- 
5937 (Gavriil et al. 2002; Gavriil et al. 2006) and in one 
out of the four bursts from XTE J1810-197(Woods et al. 
2005). We have reanalyzed all these burst spectra con- 
sistently as for 4U 0142+61. In Figure 13 we plot the 
spectra of all AXP bursts with likely emission lines in 
their spectra. Notice that all have broad features that 
occur between 13 and 14 keV. There is possible evidence 
for features at and ~21 keV in the 4U 0142+61 burst; 
there may be hints of these features in some of the other 
burst spectra as well. 

If the features are interpreted as of a cyclotron origin, 
it is puzzling that three sources with spin- down- inferred 
magnetic fields that span a range of a factor of three pro- 
duce features at similar energies. The similarity in energy 
suggests a mechanism that is not very sensitive to the ex- 
act value of the magnetic field. Moreover, it is unclear 
wdiy such features are seen only in some bursts. Thus far, 
all AXP burst spectral features have been observed with 
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Year 

2006 2007 



Modified Julian Day 

Fig. 12. Panel 1: Frequency vs. time plot of local linear 
ephemerides with uncertainties. The linear trend corresponding to 
the local ephemeris directly before the active phase is subtracted 
from all local ephemerides. Panel 2: Timing residuals after sub- 
tracting the ephemerides shown in panel 1. Panel 3: Frequency vs. 
time plot (with uncertainties) of the best-fit linear ephemerides in 
the pre-active phase region and in the post burst 2 region. The 
same trend as in panel 1 was subtracted. Panel 4: Timing resid- 
uals after subtracting the ephemerides shown in panel 3. Panel 5: 
Frequency vs. time plot (with uncertainties) of the best-fit long- 
term ephemeris which includes a glitch at MJD 53809 followed by 
a fast exponential recovery. The same trend as in panel 1 was sub- 
tracted. Panel 6: Timing residuals after subtracting the ephemeris 
shown in panel 5. 
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Fig. 13. — Burst spectra of all observed AXP bursts with signif- 
icant emission lines, as seen by RXTE. Panel 1A: Burst spectrum 
of 4U 0142+61 burst 6. The dotted line indicates the continuum 
(blackbody) component of the best-fit model. Panel IB: Residuals, 
expressed in terms of the change in y 2 , after subtracting only the 
continuum component of the best-fit model. Panels 2 A and 2B: 
Same but for burst 4 of XTE J 1810-197 (see Woods et al. 2005). 
Panels 3A and 3B: Same, but for burst 3 of IE 1048.1-5937 (see 
Gavriil et al. 2006). Panels 4A and 4B: Same, but for burst 1 of 
IE 1048.1-5937 (see Gavriil et al. 2002). 


RXTE. Their absence from many other RXTE - observed 
bursts, together with the absence of any response feature 
in the PC A near 13-14 keV, is supporting evidence for 
their veracity. Nevertheless, it is important to confirm 
these features with other instruments. Observations of 
AXPs during active phases using the Swift XRT or the 
large- area LAXPC aboard the planned Astrosat mission 
would be particularly useful in this regard. 


4.2. The Net Spin-Down Event in 4U 0142+61 

The timing glitch reported on in §3.5 had recovery 
fraction, defined as Q = Ai^/(A va ± A*/), equal to 
1.07 ±0.02. Q > 1 implies that the net frequency change 
after the transients have decayed is negative (see panel 3 
of Figure 12 and Table 4). Indeed, the net A vjv ob- 
served is — 1.1 ± 0.1 x 10”'. Such a Q > 1 has not 
previously been reported in any AXP but has been seen 
following magnetar-like radiative behavior in one high-5 
rotation-powered pulsar (Livingstone et al. 2009). Note 
that Morii et al. (2005) reported a candidate glitch in 
1999 from this source. If such a glitch occured, based 
on the ephemerides reported in Dib et al. (2007), the 
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fractional frequency change, after any recovery, would 
have been +1.9(2) xl0~ 7 — 8.6(2) xlO -7 depending on the 
glitch epoch, but manifestly positive. 

In the standard model for glitches in rotation-powered 
pulsars, the neutron-star crust contains superfluid neu- 
trons that rotate faster than the hulk of the surrounding 
matter (see, e.g., Anderson h It oh 1975; Anderson et al. 
1982; Alpar et al. 1984; Alpar & Pines 1993). This is ar- 
gued to be a result of the fact that the external magnetic 
torque acts on the crust and coupled core components 
only, with the uncoupled super fluid unaffected. The su- 
perfluid’s angular momentum resides in quantized vortex 
lines whose density is proportional to angular velocity. 
The vortex lines are suggested to be pinned to crustal nu- 
clei, and suffer strong forces due to the angular velocity 
lag that builds between crust and superfluid. A glitch in 
this picture is a sudden unpinning and outward motion of 
vortex lines, with a transfer of angular momentum from 
superfluid to crust. In a magnetar, unpinning may occur 
due to the strong internal magnetic field as it deforms 
the crust plastically or cracks it violently (Thompson & 
Duncan 1996). 

In the 4U 0142+61 event we describe in §3.5, if the 
standard glitch model is roughly correct and applies here, 
then some regions of the stellar superfluid were origi- 
nally spinning slower than the crust. Then the tran- 
sient increase in frequency would be a result of transfer 
of angular momentum first from a more rapidly rotat- 
ing region, with a subsequent angular momentum drain 
from the crust to a more slowly rotating region. As 
argued by Thompson et al. (2000), regions of slower- 
rotating superfluid can occur in magnetars, because the 
superfluid vortex motions are governed not be spin- down- 
related forces, but by advection across the stellar surface 
by the deforming crust. Those authors show that the 
number of vortex lines per unit surface area of crust 
can increase or decrease depending on the crustal mo- 
tion relative to the stellar rotation axis. They invoke 
this possibility to explain a possible “anti-glitch” seen in 
SGR 1900+14. However, that event was orders of magni- 
tude larger than what we have observed in 4U 0142+61, 
with A vjv ~ 10 -4 (Thompson et al. 2000), and also was 
likely associated with the extremely energetic flare seen 
on 1998 August 27 (Mazets et al. 1999). 

Independent evidence for slow crustal deformations in 
4U 0142+61 that could result in the regions of slower- 
rotating superfluid required to explain the spin-down 
may come from the long-term flux and pulse profile evo- 
lution previously reported for this source by Dib et al. 
(2007) and shown again here in Figures 5 and 6. Al- 
though such variations could also be magnetospheric in 
nature, possibly due to twisting of the magnetic field 
(e..g. Thompson et al. 2002; Beloborodov & Thompson 
2007), problems with this interpretation exist, as dis- 
cussed by Dib et al. (2007). For example, if magneto- 
spheric, it is puzzling that- most if not all of the changes 
in pulse profile are seen below 4 keV, with none above 
6 keV. On the other hand, slow evolution of the sur- 
face emission (e.g. Ozel & Giiver 2007), hence structure, 
is perhaps more consistent with low-energy evolution. It 
remains to be seen if detailed modeling can reproduce the 
sort of slow pulse profile changes we have observed (Fig. 
5), quite apart from the sudden appearance of new peaks, 


hence greater harmonic structure, near bursts. We note 
that the latter has also recently been seen during periods 
of activity in IE 1048.1-5937 (Tam et al. 2008). 

Previously Kaspi et al. (2000) noted that the existence 
of AXP glitches having properties comparable to those 
seen rotation-powered pulsars was not in itself evidence 
for AXPs being magnetars. This is because the stan- 
dard glitch model applies regardless of the mechanism 
by which the crust slows down; in principle a glitch can 
occur in an accretion-powered neutron star, since in this 
case too there should develop a crust/superfluid angu- 
lar velocity lag. However only in the context of magne- 
tars has the possibility of more slowly rotating superfluid 
been suggested (Thompson et al. 2000). The net spin- 
down in 4U 0142+61 may thus add to the already large 
amount of evidence against any accretion-powered origin 
for 4U 0142+61. 

4.3. The Active Phase of 4 U 0142+61 and Other AXP 
Outbursts 

In many ways, the 2006-2007 active phase of 
4U 0142+61 is similar to other phases of activity seen 
in AXPs: it was punctuated by short bursts, pulse pro- 
file changes were seen, and it was accompanied by a sig- 
nificant rotational anomaly. However the 4U 0142+61 
activity is unique in one interesting way: the pulsar did 
not suffer a large, long-lived pulse flux increase at the 
beginning of the phase. Indeed its pulsed flux (Fig. 9), 
apart from very near bursts, has remained relatively sta- 
ble. We note that observations with focusing X-ray tele- 
scopes may reveal some phase- averaged flux variations 
(see Gonzalez et al. 2009); indeed pulsed fraction has 
been shown to be inversely correlated with total flux in 
IE 1048.1-5937 (Tiengo et al. 2005; Tam et al. 2008), 
rendering pulsed flux variations smaller relative to to- 
tal variations. Still, for 4U 0142+61, such a correlation 
would have to conspire to render the pulsed flux steady. 
This seems unlikely, however only focusing telescope ob- 
servations can rule this out. 

Other AXP radiative outbursts have been suggested 
as arising from large magnetospheric twists, with associ- 
ated magnetospheric currents returning to the stellar sur- 
face and heating it, resulting in increased X-ray emission 
from the source (e.g. Thompson et al. 2002; Beloborodov 
L Thompson 2007; Beloborodov 2009). Such twists are 
thought to represent a release of magnetic energy and 
helicity from the internally wound-up magnetic field. 

On the assumption that the pulsed flux is a reasonable 
proxy for the total flux for 4U 0142+61, here, we find no 
clear evidence for significantly increased X-ray emission 
on time scales longer than a few minutes. Thus any large- 
scale magnetospheric twist scenario is problematic. As 
originally showed by Thompson et al. (2002), the X-ray 
luminosity due to returning currents reheating the sur- 
face in a significantly twisted magnetosphere in general 
are comparable to that produced from internal processes; 
this is clearly not observed in the 2006-7 active period of 
4U 0142+61. 

Instead, for this source, long-term evolution of the 
crust, driven presumably internally by field decay, and re- 
sulting in multiple unstable configurations though with- 
out any large scale magnetospheric twists, could result 
in sudden cracking and local rearrangements. This could 
be accompanied by bursts and profile changes, as well as 
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with vortex line shifting. Why such surface motion does 
not produce significant magnetospheric twists is, how- 
ever, puzzling, given that field lines are thought to be 
anchored in the crust; perhaps only for large motions do 
field lines become sufficiently twisted for enhanced X-ray 
emission to be produced. 

Dib et al. (2008) showed that other AXP glitches have 
been unaccompanied by radiative changes. For exam- 
ple, AXP IE 1841-045 has glitched multiple times yet its 
pulsed flux remains very steady. Such radiatively “silent” 
glitches may be a result of internal activity that does 
not result in significant twisting of outer magnetosphere 
field lines, where as radiatively “loud” glitches could be 
those for which the magnetosphere is impacted. We note 
in fact that thus far, the data are consistent with all 
AXP radiative outbursts being accompanied by timing 
anomalies. This may also be true of SGRs, though in 
those cases, timing anomalies are harder to establish be- 
cause of the difficulty in achieving phase- coherent timing 
in quiescence. 

5. SUMMARY AND CONCLUSIONS 

We have reported on anomalous X-ray pulsar 
4U 0142+61 entering an active phase which was preceded 
with by a long-term increase in pulsed flux. The active 
phase, which commenced in 2006 March, consisted of a 
timing anomaly that can be described as a net “anti- 
glitch,” that is, a net spin- down following an initial spin- 
up that decayed on a time scale of 17 days. Following 
the glitch, we detected six bursts from the pulsar, the 
first ever observed from this source. Despite 10 yr of 
RXTE monitoring, the bursts all occurred in the nar- 
row time span between 2006 April 6 and 2007 Febru- 
ary 7. The sixth and largest burst had a unusual spec- 


trum that cannot be described by any simple continuum 
model. Rather, it can be well fit if a prominent spec- 
tral feature at ~14 keV is included, as has been reported 
in other AXP bursts, as well as weaker features at ~8 
and ~21 keV. The pulse profile of the source changed 
from double- to triple- peaked near the bursts, and un- 
derwent considerable evolution otherwise. At the burst 
epochs the relative intensity of the three peaks signifi- 
cantly varied. The pulse profile is now relaxing to its pre- 
active phase morphology. Most aspects of 4U 0 142+6 l’s 
emission changed during the active phase, with the no- 
table exception of the pulsed flux (except near bursts). 
This argues against this event being associated with a 
sudden magnetospheric twist, as has been invoked for 
other AXP activity, and is suggestive of crustal evolu- 
tion driven internally by the large magnetic field, though 
without significant magnetospheric twisting. We suggest 
that other, radiatively silent AXP glitches have a simi- 
lar origin, whereas radiatively loud AXP timing events 
occur when the crustal motions cannot avoid significant 
twisting of the magnetic field lines. 
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